Dense stellar matter with trapped neutrinos under strong magnetic fields 



Aziz Rabhi^'^'* and Constanga Providencia^' ^ 

^ Centra de Fisica Computacional, Department of Physics, 
University of Coimbra, 3004-516 Coimbra, Portugal 
^ Laboratoire de Physique de la Matiere Condensee, 
Faculte des Sciences de Tunis, Campus Universitaire, Le Belvedere-1060, Tunisia 

(Dated: September 6, 2009) 

We investigate the effects of strong magnetic fields on the equation of state of dense stellar 
neutrino- free and neutrino-trapped matter. Relativistic nuclear models both with constant couplings 
(NLW) and with density dependent parameters (DDRH) and including hyperons are considered . 
It is shown that at low densities neutrinos are suppressed in the presence of the magnetic field. The 
magnetic field reduces the strangeness fraction of neutrino-free matter and increases the strangeness 
fraction of neutrino-trapped matter. The mass-radius relation of stars described by these equations 
of state are determined. The magnetic field makes the overall equation of state stiffer and the 
stronger the field the larger the mass of maximum mass star and the smaller the baryon density at 
the center of the star. As a consequence in the presence of strong magnetic fields the possibility 
that a protoneutron star evolves to a blackhole is smaller. 

PACS numbers: 26.60.-l-c, 12.39.Ba, 21.65.+f, 97.60. Jd 



I. INTRODUCTION 

Neutron stars with very strong magnetic fields of the order of 10^'' — 10^^ G are known as magnetars [1-3]. They 
are believed to be the sources of the intense 7 and X rays (for a review refer to [4]). Although until presently only 
16 magnetars have been identified as short 7-ray repeaters or anomalous X-ray pulsars [5], according to Ref. [6], a 
fraction as high as 10% of the neutron star population could be magnetars. Magnetars are warm, young stars, ^ 1 
kyear old. 

The study of dense stellar matter under strong magnetic fields is, therefore of great interest. Our knowledge of 
neutron star composition and structure is still uncertain [7] . For densities below twice normal nuclear matter density 
(po ~ 0.153 fm^'^), the matter consists only of nucleons and leptons. For baryon densities above 2/9o, the equation 
of state (EOS) and the composition of matter are much less certain and the strangeness degree of freedom should be 
taken into account either through the onset of hyperons, kaon condensation or a deconfinement phase transition into 
strange quark matter. The presence of hyperons in neutron stars, which has been studied by many authors [8-13], 
tends to soften the EOS at high density and lower the maximum mass of neutron stars [8-14]. Cold dense stellar 
matter is neutrino free but the protoneutron star formed after a supernova explosion, with an entropy per baryon of 
the order of 1 to 2, contains trapped neutrinos. After 10 to 20 s, the star stabilizes at practically zero temperature 
and no trapped neutrinos are left [15]. 

In this article, we focus on the properties of hyperonic matter, which is composed of a chemically equilibrated and 
charge-neutral mixture of nucleons, hyperons, and leptons. The meson-hyperon couplings play an important role in 
the determination of the EOS and the composition of hyperonic matter. In the presence of strong magnetic fields, 
the pressure and composition of matter can be affected significantly [16]. In Rcf. [16], the authors have investigated 
the effects of strong magnetic fields on the properties of neutron star matter including hyperons, but no neutrino 
trapping was considered. It was found that the EOS at high density could be significantly affected both by Landau 
quantization and by the magnetic moment interactions, but only for field strength B > 5 x 10^® G. It was also shown 
that the threshold densities of hyperons can be significantly altered by strong magnetic field. Similar conclusions 
were obtained in [17] where the strangeness was included through an antikaon condensation or in [18] where not only 
hyperons but also the strange mesons a* and (f> were included in the EOS. 

It is very interesting to investigate the influence of strong magnetic fields on hyperonic matter when neutrinos are 
trapped. Although, in reality we should consider warm matter, it was shown [19] that the effect of the inclusion of 
trapped neutrinos is more important than the temperature effect. We expect, therefore, that the main conclusions 
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taken for T = will still be valid at finite temperature. For matter without a strong magnetic field, it has been shown 
in Ref. [15, 19] that the EOS with trapped neutrinos are harder than the neutrino-free EOS, and, the maximum 
baryonic allowed mass of a stable star is higher when neutrinos are trapped. As a consequence, after the neutrino 
outflow has occurred the most massive stars may decay into a blackhole. 

The effect of the magnetic field on the structure and composition of a neutron star allowing quark-hadron phase 
transition with trapped neutrinos has been studied in Ref. [20] . They have concluded that the strong magnetic field 
makes the overall equation of state softer giving rise to much smaller maximum mass configurations which do not 
favor the formation of low mass blackholes. Their conclusion about the reduction of the maximum mass configuration 
is in contrast with the result from other works [16, 21] where an increase is predicted due to the positive contribution 
of the magnetic field pressure to the total EOS. In Ref. [20] the magnetic field pressure enters the total EOS as a 
negative contribution, making the overall EOS softer. 

This work is organized as follows: we make a brief review of the formalism used for the hadron matter with and 
without trapped neutrinos. Next we present and discuss the results obtained for the equation of state of stellar matter 
and the mass/radius properties of the corresponding compact star families, for several values of the magnetic field. 
At the end we will draw some conclusions. 



II. THE FORMALISM 



For the description of the EOS of neutron star matter, we employ a field-theoretical approach. The Lagrangian 
density of the relativistic TW model [22, 23] is given by 



nAi^** - qbliJ.A'^ -mb + TabO' 



2^ 2'^ 2'^^ 4 

1 
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- V^i^^. + Jm^^p^-ip'^^P^. (1) 



where ^Pt, and ^pi are the baryon and lepton Dirac fields, respectively. The index b runs overs the eight lightest baryons 
n, p, A, S~, E°, E+, S~ and E!° (neglecting the Q~ and the A quartet, which appear only at quite high densities, 
does not qualitatively affect our conclusions), and the sum on I is over electrons and muons (e~ and a, ui, and 
p represent the scalar, vector, and isovector-vector meson fields, which are exchanged for the description of nuclear 
interactions and A^^ = (0,0, Pa;,0) refers to a external magnetic field along the z-axis. The baryon mass and isospin 
projection are denoted by mi, and t^^ , respectively. The mesonic and electromagnetic field strength tensors are given 
by their usual expressions: fi^ji^ — d^oju — d^tOfj,, P^,y = dfj,pi, — duPp., and Ff^i, — d^A,j — d^A^. The baryon anomalous 
magnetic moments (AMM) are introduced via the coupling of the baryons to the electromagnetic field tensor with 
^Aif — t [7m '7i^] ^^'^ strength Kb = {pb/fJ'N — qbmp/mb). We neglect the AMM of the leptons in this work, because 
their effect is very small as shown in Ref. [24]. For the electromagnetic field, only frozen- field configurations are 
considered and thus there is no associated field equation. The density dependent strong interaction couplings are 
denoted by F, the electromagnetic couplings by q and the baryons, mesons and leptons masses by m. The parameters 
of the model are the nucleon mass M = 939 MeV, the masses of mesons rria,, mp and the density dependent 
coupling parameters which are adjusted in order to reproduce some of the nuclear matter bulk properties and relations 
with the Dirac-Brueckner Hartree-Fock (DBHF) calculations [25], using the following parametrisation 

Ti{p)=Ti{po)fi{x), i = a,LJ (2) 

with 

l + bi{x + dif 

fi{x)=ai , ^,.2 (3) 
l + Ci[x + di) 

where x = p/po and 



^p{p) = rp(po) exp [-ap{x - 1)] 



(4) 
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TABLE I: Static properties of the baryons considered in this study. The mass, electric charge and strange charge of the baryon 
b are denoted by mt, Qh and q^, respectively. The baryonic magnetic moment is denoted by and the anomalous magnetic 



moment by Kb 


= (Atb//ijv — qbmp/mb), where /jat 


is the nuclear ma{ 


jneton. 




UOiLyKJlL llcLllltP 


Mass 


Charge 


o 111 diii^ciicoo 


A/Tsi CTTif^'hir' mrtTTifin't' 




\) 


(McV) 


96(e) 




f^bf I^N 




p 


938.27 


1 


u 




1 7Q 


n 


939.56 








-1.91 


-1.91 


A" 


1115.7 





-1 


-0.61 


-0.61 


E+ 


1189.4 


1 


-1 


2.46 


1.67 


S° 


1192.6 





-1 


1.61 


1.61 


E" 


1197.4 


-1 


-1 


-1.16 


-038 




1314.8 





-2 


-1.25 


-1.25 




1321.3 


-1 


-2 


-0.65 


0.06 



TABLE II: Saturation properties of nuclear matter and the nucleon-meson coupling constants for the GMl and TW models. 
The relative hyperon-meson coupling constants used in the calculation are also given. 
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with the values of the parameters rrii, F;, a;, 6,;, Ci and di, i — a,U!, p given in [23]. The meson- hyperon couplings are 
assumed to be fixed fractions of the meson- nuc;leon c;onplings, = XiH^iN^ where for each meson i, the values of 
XiH are assumed equal for all hyperons H. The values of XiH are chosen to reproduce the binding energy of the A at 
nuclear saturation as suggested by Glendenning and Moszkwoski and are given in Table II. 

For GMl model, we add to the Lagrangian density, Eq. (1), the scalar meson self-interaction terms 

where gi = Ti, b and c are two dimcnsionless parameters. 

The field equations of motion follow from the Euler-Lagrange equations. From the Lagrangian density in Eq. (1), 
we obtain the following meson field equations in the mean-field approximation 



X] r^TbPb = r^AT ^ x^bpt (5) 

b b 

'^l^^ = X] ^'^bPb = ^ojN ^ X^bPb (6) 

b b 

^Ip° = X!^'''''^3i,P6 = rpAT^XpfcTs^Pb (7) 
6 6 

where a = (a), ui^ = (w^) and = (p'^) are the nonvanishing expectation values of the meson fields in uniform 
matter. 

The Dirac equations for baryons and leptons are, respectively, given by 

- qb^^A'' -ml-^o + TpT^.p'' + S^) - ^/ijvKfea^.-F''"] *6 = (8) 

{il^^di" -qa^A" -mi)^j:i = (9) 

where ml = rrib — T^a is the effective mass of baryon species b. For neutrino-free stellar matter consisting of a 
/3-equilibrium mixture of baryons and leptons, the following equilibrium conditions must be imposed: 



Mb = 96^71 - 9; Me, 



(10) 
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which are equivalent to 

= Ma = Me" = Me" , 

fJ'P = ME+ = IJ'n — IJ-ei 
Ms- = MH- = Mn + Me, 
M/:* = Me, (11) 

where /U, is the chemical potential of species i. The electric charge neutrality condition is expressed by 

6 I 

where is the number density of particle i. If trapped neutrinos are included, we replace — > Me ~ Mi^e the above 
equations, 

M6 = <lbl^n- qi {iJ-e - Mi^e ) • (13) 

For leptons we have 

Mm-Mi/^ =Me-M!^e- (14) 

In the above equations /i^, are the chemical potential, respectively, of the muon and electron neutrinos. The 
introduction of additional variables, the neutrino chemical potentials, requires additional constraints, which we supply 
by fixing the lepton fraction Y^e = + = 0.4 [15, 26]. Also, because no muons are present before and during the 
supernova explosion, the constraint Yl^ = 1^ + yj^^ = must be imposed. 

The energy spectra for charged baryons, neutral baryons and leptons (electrons and muons) are given by 



^1,8 



El. 



= Jk^,+ (^^mf + 2u\qb\B-s,iNKbB^ + r<^6W° + T3,rp6M° + (15) 



= ^k^,+mf + 2u\qi\B (17) 



where u = n + ^ — sgn{q)^ = 0, 1, 2, . . . enumerates the Landau levels (LL) of the fermions with electric charge q, the 
quantum number s is +1 for spin up and —1 for spin down cases, and the rearrangement term is given by 



V dp dp dp 

1 dT^N 2 2, 1 dTpN 22 1 dr^N 22 /1 o\ 

= — a^r^^'^o + f — a^r^ppo - f — ^"^-^^ • (18) 

i-uiN Cfp I pN dp i crJV Op 

For the charged baryons, we introduce the efFcctivc mass under the effect of a magnetic field 

ml = ^Jmf + 2v\qh\B - sp,NKbB, (19) 
the expressions of the scalar and vector densities are, respectively, given by [27] 

27r2 A. ^rnl- + 2i.\qb\B ml 

= 1^ E E'^k.s- (20) 

v=0 s 

where kp^^ is the Fermi momenta of charged baryons b with quantum numbers v and s. The Fermi energies Ep are 
related to the Fermi momenta kp^^ by 

{k''F,.,sf = (e'f)' - {<f ■ (21) 
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For neutral baryons the Fermi momenta is denoted by fc^^, and the Fermi energy Ep is given by 



{k 



F,s) - {E%) 



'6 J 



with 



The scalar and vector densities of the neutral baryon b are, respectively, given by 
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(23) 
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The vector densities for leptons are given by 
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where fc^ ^ is the lepton Fermi momenta, which are related to the Fermi energy Ep by 

ikKu,sf ^ (Epf - ml l = e,i2, 



(24) 



(25) 



(26) 



with fhf = mf + 2i'\qi\B. The summation in i/ in the above expressions terminates at fmax, the largest value of v for 
which the square of Fermi momenta of the particle is still positive and which corresponds to the closest integer from 
below. For neutrino- trapped matter, the neutrino density is given by 



fhp 

6^' 



The chemical potentials of baryons and leptons are defined as 

fxb = £^ + r^ftwO + rp6T3,p° + 



(27) 

(28) 
(29) 



We solve the coupled Eqs. (5)-(12) self-consistently at a given baryon density p = ^i,Pb presence of strong 

magnetic fields. The energy density of neutron star matter is given by 



2'^Ipo, 

b l=e,fj, 

where the energy densities of charged baryons and leptons have the following forms 
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while those of uncharged baryons are given by 
1 
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The pressure of neutron star matter can be obtained by 

-Pm = ^ IJ-iPi - £m = IJ-n^ Pb - (33) 
i b 

where the charge neutrality and /3-equihbrium conditions are used to get the last equality. If the stellar matter 
contains neutrinos trapped, their energy density and pressure contributions, respectively. 



87r2' 

should be added to the stellar matter energy and pressure. The total energy density and pressure of the system 
includes the contribution of the magnetic field, 

^2 

P = Pm + ^. (35) 

With the obtained EOS, the mass-radius relation and other relevant quantities of neutron star can be derived by 
solving the Tolman-Oppenheimer-Volkoff (TOV) equations. 



III. RESULTS AND DISCUSSION 



A. Uniform magnetic field 

In this section we consider that the external magnetic field is constant. The magnetic field will be defined in units of 
the critical field = 4.414 x lO""-^ G, so that B = B* B^. In order to study the effect of strong magnetic fields on the 
structure of hyperonic matter we use two different relativistic mean-field approaches: the GMl parametrisation of the 
NLW models [28] , and the TW parametrisation of the density-dependent relativistic hadronic (DDRH) models [22, 23] , 
given in Table II. We include the baryonic octet in the EOS and choose two sets of hyperon-meson coupling constants 
given in Table II. The static properties of the baryons considered are listed in Table I. 

In Fig. 1, the EoS obtained with both relativistic mean- field models GMl and TW are displayed for B* = 
0,10^, and 3 x 10^, for the relative hyperon coupling constant x„ = 0.6 (thick lines) and x^r = 0.8 (thin lines). 
In Fig. 1 (a) and Fig. 1 (c) the EOS does not include the contribution form the AMM and in Fig. 1 (b) and Fig. 1 
(d) the AMM was included. The kink on each one of the curves identifies the onset of hyperons. Both for TW, which 
gives a softer EOS, and for GMl, the EOS at higher densities arc harder for x„ = 0.8 than for x„ = 0.6 because a 
larger u oi p coupling constant makes the onset of hyperons only possible at larger densities. This is also valid in the 
presence of a strong magnetic field. The effects of the AMM are only noticeable for the stronger magnetic field. It 
is seen that, in the absence of the AMM and in presence of a strong magnetic field, the EOS becomes softer for the 
smaller densities and harder at the larger densities due to the Landau quantization which affects charged particles. 

The effect of magnetic field on the EOS, when trapped neutrinos are considered, is seen in Fig. 2 for GMl and TW 
models, with x^ = 0.6 and without/with AMM. For B = the presence of neutrinos makes the EOS harder [15, 19]. 
A strong magnetic field without including the AMM, makes the EOS harder at high densities if neutrino trapping is 
enforced. Moreover, if the AMM are included the EOS become even harder. At low densities, however, the magnetic 
field softens the EOS even when the AMM are taken into account. We remark that for the largest field considered and 
taking into account AMM, it is not possible to get an EOS at low densities for the lepton fraction Yl = 0.4. In Fig. 3 
we plot the baryon density threshold, at which the neutrino chemical potential vanishes, as function of the magnetic 
field, for GMl and TW models using x^ = 0.6 and without /with including AMM. Below this threshold density, it is 
not possible to impose the lepton fraction Yl = 0.4. The threshold density increases very fast with the increase of the 
magnetic field, e.g. in Fig. 3(a) for B* = 5 x 10^ corresponds a threshold density equal to 0.5po for GMl and l.'Sbpo 
for TW. The kink at p ~ 2.5po corresponds to the hyperon onset. However, at the surface the measured magnetic 
fields of magnetars are at most B* ^ 10^ and it is probable that even if the magnetic field is stronger in the interior 
it will be weaker than the threshold values given in Fig. 3. 

In Fig. 4 we show, for several values of magnetic field without including the AMM, for GMl and TW models and 
for the two values of the hyperon-meson coupling constants, the strangeness fraction defined as 

HP _ Efc \Qs\ Pb 
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where 5^ is the strange charge of baryon 6, and is listed in Table I. The strangeness onset occurs around 2po and it 
has almost reached 0.275 (0.24) for GMl and 0.23 (0.18) for TW, neutrino-free matter and x„ = 0.6 (xa = 0.8). For 
5 = and for the neutrino-free matter, increasing the hypcron-meson couplings decreases the strangeness contents, 
and this suppression of strangeness is stronger in TW than GMl. However, if the neutrino trapping is imposed the 
strangeness fraction is smaller and a stronger reduction occurs for GMl and for Xc — 0.8. 

For neutrino-free matter, the magnetic field reduces the strangeness fractions below a critical density that can be 
quite high, 6po or larger for the fields shown. This effect starts to be detected already for B* = 10^, which represents 
the threshold for the magnetic field effects to become significant, independently of the EOS [16]. However, the magnetic 
field enhances the strangeness fraction of neutrino-trapped matter. If the AMM are included the strangeness fraction 
behaves in the same way. 

In Figs. 5, 6, 7, and 8, we present the particle fractions Yi — p\l p as function of the baryon density p/ po^ obtained 
with TW model using two values of hyperon-meson coupling constants and for several values of the magnetic fields. 
Figs. 5 and 7 are for neutrino-free matter and Figs. 6 and 8 for neutrino-trapped matter. For GMl we get similar 
results. 

First, we study the results obtained with the hyperon-meson coupling constant equal to = 0.6. For B — Q 
neutrino free matter, shown in Fig. 5 (a), the proton fraction rises quickly with increasing density and reaches ~ 0.1 
around 2po before the appearance of hyperons. is the first hyperon to appear. In Fig. 5 (b) and (d), we show 
the results for B* = 10^ and B* — 2> x 10^. At low densities, the fraction of nucleons and leptons are significantly 
affected by the magnetic field. The Landau quantization increases the proton abundance, and, therefore, the electron 
abundance due to the charge neutrality. In Fig. 9 we show for TW with x^, ~ 0.6 the neutron, proton, electron and 
neutrino chemical potentials as a function of density for different values of the magnetic field. It is seen that the proton 
and neutron chemical potentials decrease with the increase of B. This is also the general tendency of leptons except 
for the larger densities. In the absence of AMM the reduction of the neutron chemical potential with an increase 
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FIG. 2; EOS for TW and GMl models and for Xa = 0.6, without (a) and (c) and with (b) and (d) AMM. Thick hnes for 
neutrino-trapped matter and thin hnes for neutrino-free matter. 



of the magnetic field is due to a reduction of the isospin asymmetry. The AMM will further reduce the neutron 
chemical potential. Landau quantization is the main explanation for a reduction of the proton chemical potential. As 
a consequence the threshold for the appearance of hyperons like Tr and A will occur at larger densities [16]. However, 
at densities above ~ Spo the threshold densities for S°, E+, S~ occur at smaller densities probably due to the smaller 
effective masses of these hyperons at large densities with strong magnetic fields. The inclusion of AMM, Fig. 5 (c) 
and (e), produces an even larger hyperon suppression at the lower densities. 

We consider now matter with trapped neutrinos. For _B = 0, Fig. 6 (a), the proton and electron fractions are ^ 0.3 
from low densities until the appearance of hyperons at ^ 3po and, therefore, the neutrino fraction is ^ 0.1 since 
we are imposing 1^ = 0.4. A is the first hyperon to appear for B=0. For B* = 10^ (see Fig. 6 (b)), the proton 
fraction is larger at the smaller densities, ~ 0.4, decreases to 0.3 for p ~ po and stays approximately constant until 
the appearance of hyperons at ~ 3po- The neutrino fraction rises quickly with increasing density and reaches ^0.1 
around po and keeps, roughly, this value, ~ 0.1. A is still the first hyperon to appear in this case. However, for 
= 3 X 10^ (see Fig. 6 (d)), the proton fraction decreases slightly starting at ~ 0.4 at lower densities until the 
appearance of hyperons at ^ 3po- is the first hyperon to appear almost at the same density as A. The neutrino 
fraction is almost zero at subsaturation densities, and rises, quickly, with increasing density, to ~ 0.2 around 6po- The 
net effect of magnetic field is a neutrino suppression due to the larger proton and therefore electron fractions and to 
move the threshold density of the negatively charged baryons (e.^. E~) to larger densities and the positively charged 
baryons to lower densities. This is due to a decrease (increase) of the chemical potentials of negatively charged baryons 
(positively charged baryons) when the neutrino chemical potential is taken into account. In particular the onset of 
the S~ occurs for densities larger than 8po- The neutral baryons are not affected by the presence of neutrinos. In the 
presence of a strong magnetic field the overall effect of the neutrinos is smaller than for B=0. 

For x„ = 0.8, we obtain similar results and the main differences are: at high densities the onset of S~ occurs at 
p < 8po, the onset of the S~ and A occurs at larger densities and hyperon fractions are smaller. 

The effect of the magnetic field on the onset of hyperons, in both cases for neutrino free and neutrino trapped 
matter, is clearly shown in Tables HI and IV, where we give the threshold density of the first hyperon, for TW and 
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GMl models and for Xa = 0.6, and 0.8, without including the AMM (Table III)and including the AMM in (Table IV). 
The main conclusions for neutrino-free matter are: a) the hyperon onset occurs with the appearance of the E~ meson, 
both with and without AMM; b) the onset density generally increases with the increase of the magnetic field but for 
GMl it may decrease for B* — 10^; c) the inclusion of AMM reduces the onset density of hyperons; d) the hyperon 
onset occurs at larger densities for x„ = 0.8. For neutrino-trapped matter we conclude that: a) the onset density 
of hyperons is generally get smaller when the magnetic field magnitude increases; b) for small fields A is the first 
hyperon to appear but for large fields it is either S~ or S~ and A together. 

In order to better understand the effect of the magnetic field on the neutrino trapping, we show in Fig. 10 the 
fraction of neutrinos for several values of magnetic field. For i? = 0, the neutrino fraction decreases at low densities 
and starts to increase after the onset of hyperons because with the hyperon onset the electron fraction gets smaller. 
For B* = 5 X 10"* and B* = 10^, the neutrino fraction is affected by the Landau quantization of electrons, and 
oscillates around the B = results, when AMM are not included. At low densities the main effect of magnetic field is 
the suppression of neutrinos. For B* > 10^, the neutrino fraction is lowered at low densities and enhanced at higher 
densities. However, when the AMM are included, the abundance of neutrinos is only slightly reduced for B* = 5x 10^ 
and B* = 10^, while for B* > 10^ the neutrino suppression is strong, e.g. for B* = 3 x 10^ there are no neutrinos 
at a density below the saturation density. We conclude, therefore, that the abundance of neutrinos in the presence of 
strong magnetic fields has a strong suppression at low densities. 

In order to study the effect of the magnetic field on the properties of the stars described by the EOS discussed 
above with neutrino trapping, we must consider a magnetic field depending on the baryon density, which is not larger 
than lO^^G at the surface. We will do this in the next section. 
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B. Baryon density-dependent magnetic field 



Since, to date, there is no information available on the interior magnetic field of the star, we will assume that the 
magnetic field is baryon density-dependent as suggested by Ref. [29]. The variation of the magnetic field B with the 
baryon density p from the center to the surface of a star is parametrized [29, 30] by the following form 



1 — exp < —(3 



Po 



(36) 



where po is the saturation density, iJ^urf jg magnetic field at the surface taken equal to lO^^G, in accordance with 
the values inferred from observations and Bq represents the magnetic field at large densities. The parameters /3 and 
7 may be chosen in such way that the field decreases fast or slowly with the density from the center to the surface. In 
this work, we will use one set of value (/3 — 0.05 and 7 = 2) allowing a slowly varying field. The magnetic field will be 
given in units of the critical field B^ = 4.414 x lO^'^ G, so that Bq = Bq B^. We further take Bq as a free parameter 
to check the effect of different fields. A detailed discussion about the variation of the magnetic field with the baryon 
density may be found in Ref. [31]. 

In this section, we will only present the results for the density-dependent relativistic model TW and we will take for 
the hyperon-meson coupling constants Xa- = 0.6. Furthermore, in all the figures, we will only show the results obtained 
without including the baryonic AMM, because for the intensity of the magnetic fields considered its contribution to 
the EOS is negligible. 

In Fig. 11, we show the total pressure P, see Eq. (35), as a function of the total energy density e for the magnetic 
field strengths B* = 0, 10^, and 3 x 10^. The results with neutrino-trapped and neutrino-free matter are plotted with 
thick and thin lines, respectively. The EOS with neutrino trapping is stiffer than the neutrino free EOS, both with 
and without magnetic field. However, when neutrino trapping is imposed the magnetic field makes the EOS softer at 
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FIG. 6: Baryon fraction versus baryon density obtained with TW model for Xa = 0.6, and neutrino-trapped matter. 



low density. 

We have also studied the baryonic and leptonic composition of the stars. In Fig. 12 we plot the lepton and the 
baryon fractions for neutrino free matter (left column) and for neutrino trapped matter (right column) for TW model. 
We show the particle fractions for the smallest and the largest magnetic fields considered in the present work. The 
vertical lines in the figures represent the central density of the star configuration with maximum mass when it lies 
within the range of densities shown. We conclude that the main effect of magnetic fields with the intensity considered 
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is the appearance at high densities of a larger leptonic fraction. Due to the contribution of the magnetic field energy 
and pressure to the total EOS of stellar matter the contribution of matter becomes smaller when the magnetic field 
increases and, therefore, the fraction of strangeness in the star is strongly reduced because the larger contributions of 
hyperons for the star come from the larger densities. 

Hadron star profiles are obtained from the EOS studied, for severals values of magnetic field, by solving the Tolman- 
Oppenheimer-Volkoff equations, resulting from Einstein's general relativity equations for spherically symmetric and 
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static stars. Although an approximation due to the presence of the magnetic field, we assume spherical symmetric 
stars. In Table V we show the values obtained for the maximum gravitational and baryonic masses of the stars, their 
radius, their central energy and baryon densities and the magnetic field at the center. The results are shown for 
the stars with and without trapped neutrinos. The mass/radius curves for the families of stars corresponding to the 
maximum mass configurations given in Table V are plotted in Fig. 13; the thick lines correspond to the stars with 
the trapped neutrinos. From Table V, we conclude that maximum mass and radius of a neutrino free star is generally 
smaller than the corresponding mass and radius of a star with trapped neutrinos. The presence of neutrinos makes 
the EOS stiffer. This difference, however, is reduced with the increase of the magnetic field. Quantitatively, in the 
field free case the difference is around 1.75 km and it decreases to 1.22 km for Bq = 5 x 10'', to 0.85 km for Bq — 10^, 
and to 0.09km Bq — 3 x 10^. This is due to the fact that the contribution of matter becomes smaller and smaller as 
the magnetic field increases, and therefore the mass and radius are not so sensitive to the specific properties of the 
stellar matter contribution to the total EOS. It has already been shown in several works [16, 21, 31] that the stronger 
the magnetic field the stiffer the EOS and therefore the larger the star radius and mass. 

An analysis of the Table V allows us to draw some interesting conclusions: a) although the presence of neutrinos 
makes the maximum gravitational and baryon masses larger, the difference with respect to the maximum gravitational 
and baryon masses of a neutrino free star decreases with the increase of the magnetic field intensity. For the largest 
field considered the star with neutrinos has a smaller mass; b) as a consequence in the presence of a very strong 
magnetic field a star will never decay into a blackhole after the outflow of neutrinos. This, however, may occur if 
no magnetic field/a weak magnetic field exists because the baryonic mass of the maximum mass star configuration 
is larger if neutrinos are trapped in the star; c) the most massive neutron star reported till recently was the PSR 
J1748-2021B with a mass 2.74±0.21 [32], which, however, was still not confirmed by other measurements. This kind 
of measurements may impose a constraint on the maximum acceptable magnetic field. For Bq = 3 x 10'^ we get at 
the center of the star B ~ 4 x lO^^G and a mass larger than 3.2 Mq. 
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IV. CONCLUSIONS AND OUTLOOKS 



In the present work we have studied the effect of a strong magnetic field on the EOS and the star properties 
when neutrinos are trapped in the star. We have used two different relativistic mean-field models, one with constant 
coupling parameters, GMl, and the other with density dependent coupling parameters, TW. We have also considered 
two sets of hyperon-coupling parameters. The main conclusions of the work do not depend either on the model or on 
the strength of the hyperon-couplings. 

The phase of trapped neutrinos in the life of a proto-neutron star occurs while the stellar matter is still warm 
[15, 26], and, therefore, a finite temperature calculation should have been done. However, we do not expect that 
temperature will change the main conclusions of the present work. In fact, in several works [15, 33, 34] it has been 
shown that the star properties such as mass and radius do not depend much on temperature. 

We have shown that a strong magnetic field suppresses the presence of neutrinos at low densities. It was also 
shown that although strangeness is suppressed by the presence of neutrinos, if the star has a strong magnetic field 
this suppression is smaller. The magnetic field affects in a different way the charged and neutral baryons and it may 
affect the order at which density they appear. For neutrino free matter, E~ is the first hyperon to appear at the 
smallest densities. However, for neutrino trapped matter, A is the first hyperon to appear except for the largest field 
considered when we may have or and A. 

We have studied the properties of stars with trapped neutrinos and strong magnetic fields: it was shown that the 
magnetic field increases the mass and radius of the most massive star configuration, and in general the radius of all 
stars with a mass larger than 0.5 Mq. The mass of observed neutron stars may set an upper limit on the possible 
magnetic field acceptable in the interior of a star. The contribution of the magnetic field to the total EOS of the star 
reduces the relative importance of the stellar matter term: the central baryonic density decreases as the field increases. 
One of the main consequences is the reduction of the strangeness fraction in the star or other exotic components such 
as a quark core [31], or kaon condensation [17]. It was also shown that the magnetic field reduces the possibility of 
formation of a blackhole after the outflow of neutrinos. 
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FIG. 11: EOS for stellar matter for TW model without/with neutrino trapping, for ~ 0.6, for several values of magnetic 
field and without AMM, using a slowly varying parametrisation of B. Thin lines correspond to neutrino- free matter and thick 
lines to neutrino-trapped matter. 
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FIG. 13: Mass-radius curve of neutron stars for several values of the magnetic field, using a slowly varying parametrisation of 
B, for TW model with/without neutrino trapping, and for = 0.6. Thin lines correspond to neutrino free matter and thick 
lines to neutrino trapped matter. 



TABLE III: The threshold density of hyperons in fm""^ when AMM is not included, for several values of magnetic field, for 
both GMl and TW models and for neutrino free and neutrino-trapped matter. We also identify, between brackets, the hyperon 
that appears at smaller densities. 



B* 


Models 




Xct = 0.6 




x„ = 0.8 


VLe =0 


Fie = 0.4 


KLe =0 


Yl. = 0.4 





TW 


0.275 (E-) 


0.441 (A) 


0.301 (E-) 


0.514 (A) 




GMl 


0.280 (E-) 


0.419 (A) 


0.327 (E-) 


0.501 (A) 


10® 


TW 


0.311 (E-) 


0.438 (A) 


0.338 (E-) 


0.516 (A, E") 




GMl 


0.301 (E") 


0.418 (A) 


0.340 (E") 


0.503 (A) 


3 X lO'^ 


TW 


0.358 (E") 


0.445 (E") 


0.385 (E") 


0.506 (A, E") 




GMl 


0.373 (E") 


0.427 (A, E~) 


0.433 (E~) 


0.496 (A) 
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TABLE IV: The threshold density of hyperons in fm~^ when AMM is included, for several values of magnetic field, for both 
GMl and TW models and for neutrino free and neutrino-trapped matter. We identify, between brackets, the hyperon that 
appears at smaller densities. 



B* 


Models 




Xa = 0.6 


Xa = 0.8 




y^e = 


YLe = 0.4 




YLe = 0.4 





TW 


0.275 (E-) 


0.441 (A) 


0.301 (E-) 


0.514 (A) 




GMl 


0.280 (S-) 


0.419 (A) 


0.328 (E-) 


0.501 (A) 


10^ 


TW 


0.291 


0.433 (A) 


0.315 (E-) 


0.500 (A) 




GMl 


0.287 (E-) 


0.409 (A, T.-) 


0.325 (E^) 


0.488 (A) 


3 X 10^ 


TW 


0.344 (E-) 


0.389 (T.-) 


0.373 (E-) 


0.457 (E-) 




GMl 


0.359 (S-) 


0.375 (E-) 


0.415 (E-) 


0.468 (E-) 



TABLE V: Compact star properties using TW model, for several values of a density-dependent magnetic field. Mmax, M^^xj 
R, Eo, p'^, and B* are the gravitational and baryonic masses, the radius, the central energy and baryon density, and the value 
of the magnetic field at the center, respectively. 







BS 


MmaxlMe] 




R [km] 


i?o[fm-*] 


u" = p^/po 


B*e 


YLe 


= 0.4 


B* = 


1.77 


1.95 


11.75 


6.35 


6.627 







= 


B* =0 


1.51 


1.70 


10.00 


7.74 


8.460 





YLe 


= 0.4 


5 X 10"* 


1.97 


2.17 


11.81 


6.10 


6.041 


4.196x10" 




= 


5 X 10"* 


1.78 


2.00 


10.59 


7.25 


7.514 


4.705x10" 


YLe 


= 0.4 


10^ 


2.31 


2.51 


12.54 


5.34 


4.740 


6.750x10* 


Y., 


= 


lO'^ 


2.21 


2.45 


11.69 


5.82 


5.303 


7.551x10* 


YLe 


= 0.4 


3 X 10^ 


3.18 


3.28 


15.76 


3.22 


2.365 


7.323x10* 




= 


3 X 10® 


3.21 


3.43 


15.67 


3.15 


2.376 


7.379x10* 



